Spin-orbit induced hole spin relaxation in InAs and GaAs quantum dots 
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We study the effect of valence band spin-orbit interactions on the acoustic phonon assisted spin 
relaxation of holes confined in quantum dots. Heavy hole-light hole (hh-lh) mixing and all the spin- 
orbit terms arising from zinc-blende bulk inversion asymmetry (BIA) are considered on equal footing 
in a fully 3D Hamiltonian. We show that hh-lh mixing and BIA have comparable contributions to 
the hole spin relaxation in self-assembled QDs, but BIA becomes dominant in gated QDs. The 
dependence of the hole spin relaxation on the QD geometry and spin splitting energy is drastically 
different from that of electrons, with a non-monotonic behavior which results from the interplay 
between SOI terms. Our results reconcile contradictory predictions of previous theoretical works 
and are consistent with experiments. 

PACS numbers: 73.21.La,71.70.Ej,72.25Rb,63.20.kd 



I. INTRODUCTION 

In the last years, the spin of holes confined in III-V 
semiconductor quantum dots (QDs) has emerged as a 
promising building block for spintronic and spin-based 
quantum information devices— As compared to electrons, 
the p-like nature of the hole orbitals leads to weaker hy- 
pcrfine interaction with the lattice nuclei, resulting in 
coherence times which hold promise for applications^ — 
As a matter of fact, demonstrations of hole spin ma- 
nipulation in QDs have been recently reported^— and 
theoretical proposals of control mechanisms are being 
proposed— i - — In this context, the study of hole spin re- 
laxation has become a subject of interest. Hole spin re- 
laxation is also important for optical applications because 
it is believed to rule the exciton spin dynamics in both 
dark-to-bright exciton transitions 1 ^— -which affect exci- 
ton storage tiniest— - and transitions within the bright 
doublet^ -which affect light depolarization^!-. 

Experimental observations in self- assembled InAs QDs 
point at hole spin lifetimes ranging from ~ 10 ps 
to Ti ~ 1 ms^— The large dispersion is partly at- 
tributed to the different relaxation mechanism involved 
in different studies. When the energy splitting between 
orthogonal spin states is small, hyperfine interaction is 
the dominant relaxation channeL 16 ' 27 In this case, the 
lifetime is strongly dependent on the degree of hh-lh mix- 
ing. If the hole state is a pure hh, as in the ground state 
of flat (quasi-2D) QDs, the hyperfine interaction takes 
an Ising form and spin relaxation is slow, but it rapidly 
increases in non-flat QDs due to hh-lh mixing. 1 ' 7 On the 
other hand, when the energy splitting exceeds the nuclear 
magnetic field, the valence band spin-orbit interaction 
(SOI) takes over as the main source of relaxation ! 16 i 27 
Long hole spin lifetimes have then been observed, reach- 
ing up to T/ 1 ~ 0.25 ms, which is only 5-10 times shorter 
than electron spin lifetimes, T-f— £ This result is encour- 
aging for the use of holes in quantum information and 
optical applications, but it is surprising because the va- 
lence band SOI interaction is known to be much stronger 
than that of the conduction bandit 



The above paradox has prompted a number of theoret- 
ical works trying to understand which factors determine 
the relaxation dynamics of single holes under magnetic 
fields 2 ^— and that of holes forming exciton a 15 i 22 ' 34 i 35 in 
quasi-2D InAs/GaAs QDs. For the relaxation to take 
place one needs a source of energy relaxation, which 
in these systems is provided by the acoustic phonon 
bath ; 23 ' 24 plus a source of spin admixture. Woods et 
al— 2 and Lu et al— £ proposed that the latter is the 
coupling between hh and Ih subbands. Other authors 
have suggested instead that the splitting between hh 
and Ih subbands in flat QDs is large owing to confine- 
ment and strain, so that spin admixture must be due 
to other SOI mechanisms. It was then proposed that 
hole SOI should have a similar origin to that of conduc- 
tion electrons, namely the BIA of zinc-blende crystals, 
which gives rise to Dresselhaus SOI terms. 2 - Bulaev and 
Loss assumed the cubic-in-fc Dresselhaus term is domi- 
nant and showed that T± could then become compara- 
ble to Tf in flat QDs— - Other studies followed this as- 
sumption and succeeded in explaining some experimen- 
tal observations ! 15 ' 26 ' 32 By contrast, Tsitsishvili et al. 
suggested that if the lateral confinement is weak, it is 
the linear-in-A; term that dominates the mixing— 5 - Last, 
Roszak et al. 34 suggested that for holes forming excitons 
it is the electron- hole (e-h) exchange interaction together 
with the strain that gives rise to hole spin admixture. 

It is worth noting that all the previous works assumed 
a dominating SOI mechanism without actually compar- 
ing with others. In addition, simplified models disre- 
garding hh-lh mixing become highly parametric, and dif- 
ferent parameters were needed to explain different ex- 
perimental observations even in the same systemi 15 ' 26 i 32 
The lack of a comprehensive study translates into many 
open questions which show that hole spin relaxation in 
QDs is still not fully understood. To name a few: (i) 
while Ref. [29| pred icts that T/ 1 decreases with the QD 
diameter, Ref. |30l predicts exactly the opposite behav- 
ior; (ii) Ref. Sjpredicted T% > T{ in the limit of 2D 
QDs, but experiments on self- assembled QDs have only 
shown T[/T[ l = 5 — 10 ~£ so that one wonders if any 
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realistic QD structure would actually permit holes relax- 
ing slower than electrons; (iii) in excitons, the role of 
e-h exchange energy is not clear: while experiments with 
self- assembled QDs have shown neglegible dependence of 
T^r^ a strong dependence has been found in colloidal 
quantum rods.— 

In this work, we aim at covering the existing gap in 
the understanding of hole spin relaxation in QDs. We 
study the hole spin dynamics considering simultaneously 
hh-lh mixing and all the different Dresselhaus SOI terms 
arising from the BIA of zinc-blende crystals, along with 
the hole-acoustic phonon coupling. All the terms are de- 
scribed within a four-band k-p formalism and 3D Hamil- 
tonians, which allows us to provide a general overview 
on the effect of the QD size and geometry dependence 
while relying on well-known bulk parameters only. In 
this way, we are able to establish the regime of valid- 
ity of previous studies which assumed a single dominat- 
ing SOI mechanism. Furthermore, we explore spheroidal 
QDs beyond the usual quasi-2D limit, thus providing the- 
oretical assessment for developing experimental research 
with spherical and prolate QDs i 18 i 36 

We find that hh-lh is the main SOI channel in prolated 
or spherical QDs, but Dresselhaus SOI has a compara- 
ble contribution in oblatcd QDs (such as self-assembled 
dots), and it becomes dominant in quasi-2D QDs with 
very weak confinement (as in electrostatically confined 
dots). The competition between SOI coupling terms and 
the energy splitting between hh and Ih leads to a non- 
monotonic dependence of T/ 1 with the QD geometry, in 
sharp contrast with the well-known case of electrons. 
This explains the opposite trends reported by different 
theoretical studies in the literature. The dependence of 
Ti on the electron-hole exchange energy we predict is 
consistent with experiments on colloidal nanorods^, but 
it suggests that two-phonon processes are relevant in self- 
assembled QDs. In prolate QDs, where the ground state 
is formed by Ih, the spin relaxation is shown to take place 
in similar timescales as for transitions between hh states. 
However, the coupling to acoustic phonons is different, 
with deformation potential interaction being the main 
mechanism even for vanishing spin splitting energy. 
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FIG. 1: Geometry of QDs with varying lateral (a) and vertical 
(b) confinement frequency. 



A. Hamiltonian 

The system Hamiltonian reads: 

H = Hfr+ H p h + Hh-ph- (f ) 

In Eq. ([TJ), Hh is the hole Hamiltonian, 

H h = H L + H BIA + V QD 1 + H Z , (2) 

where Hl is the 4-band Luttinger Hamiltonian describing 
the coupled hh-lh bands.— It includes quadratic terms in 
k only: 



Hl = —[(71 + ^72)^ - !2{klJ 2 x + klJl + h'J 



m " 2 ' ' 2 

2~f 3 ({h X ,ky}{J x , Jy} + {ky, k Z }{jy, J Z } 

{k Zl k x }{J z ,J x })] 



z^z) 



(3) 



where mo is the free electron mass, 7^ are the Luttinger 
parameters, kj = —ihdj the j component of the linear 
momentum, {A, B} = \{AB + BA) and Ji is the i-th 
component of the angular momentum corresponding to 
the quantum number J = 3/2. To get the matrix repre- 
sentation of this hamiltonian we multiply the first term 
of equation ^ by the 4x4 unit matrix and employ the 
standard matrix representation of the J — 3/2 compo- 
nents of the angular momentum. 38 We finally get: 



Hr 



fP + Q -S R \ 

-£t p-Q R 

i?t P-Q S 

\ i?t P + Q J 



(4) 



II. THEORETICAL MODEL 
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We study the spin relaxation of holes confined in zinc 
blende QDs grown along the [001] direction. The hole 
spin states are considered split energetically, for example 
by the e-h exchange interaction in excitons or any other 
source that can be viewed as an effective axial magnetic 
field. Thus, similar results can be expected for transi- 
tions between Zeeman sublevels under moderate external 
magnetic fields. 
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V3j 2 (kl-kl) + 2iV3^ 3 k x k y j (7) 
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S = —— 2v^73 (k x - ik y )k z . 
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Hbia includes the linear and the Dresselhaus SOI third 
order in k terms sSS 
2 

Hbia = -7= C k [k x {J x , J 2 - J 2 } + cp] 



6 4 i ({k x ,ky 



+ &42 ({fc a 

+ 651 ({h 



k z } J x 



>ky k z } J x 



cp) 
cp) 



k 2 }{J x ,J 2 - J z 2 } + cp) 



+ 652 (^{J„ J 2 - J z 2 } + cp), 



where C'k, 641, 642, &51 and 652 are material dependent 
coefficients and cp stands for cyclic permutations of the 
preceding terms. The matrix form of the Hamiltonian 
terms above is given in the Appendix. One can notice 
that all BIA terms provide direct mixing between hh spin 
up and spin down (J z = +3/2 and J z = —3/2) compo- 
nents except for &41, which requires the participation of 
the Ih (J z — +1/2 and J z = —1/2) components. Rashba 
SOI is neglected in this study because for holes in QDs 
it is inefficient. 31 Vqd describes the confining potential 
of the QD. We model QDs with parabolic confinement in 
the x, y and z directions, 

Vqd = -\ Xx (x 2 + y 2 ) - \ XzZ 2 (10) 

where xx an d Xz are the force constants perpendicu- 
lar and parallel to the growth direction, respectively. 
Eq. (fpQ|) allows us to simulate 3D spheroidal QDs with 
different aspect ratios, from flat (quasi-2D) to spherical 
or elongated (quasi-lD) structures, see Fig. Q] Hz is 
the Hamiltonian modeling the splitting of the hole states 
by an effective axial magnetic field, three times larger 
for heavy (\J Z \ =3/2) than for light holes (\J Z \ = 1/2). 
This field could originate e.g. from the e-h exchange 
interaction 3 ^ or a spin Zeeman effect. Then, we assume 
that: 
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(11) 



To calculate the hole states from Hh, we note that the 
diagonal terms correspond to harmonic oscillator Hamil- 
tonians: 



P + Q + V QD 



T hh,± ~hx± (z 2 + y 2 ) 

Thh,z 2 z 1 



(12) 



and 



P-Q + Vqd = T lh ^-- x ±(x 2 +y 2 ) 



1 



XzZ 



(13) 



- 2^ k j ' with i = ( hh > i = z )> fcj 



where 

(kl + k 2 ) 1 / 2 , mf = mo/(7l+72), =m /( 7l -2 7 2), 



mf = mo/ (71 - 72), and m l z = mo/ (71 + 2 72 ). This 
suggests rewriting all derivatives and coordinates of Hh 
in terms of harmonic oscillator ladder operators and then 
projecting it onto a basis formed by oscillator eigenstates. 
The problem is that Eq. (fT2j) has hh masses while Eq. (|T3]) 
has Ih masses, and hence they have different oscillator fre- 
(9) quency, oj* = [xj /"jj) 1 ' 2 ■ Because the off-diagonal terms 
of Hl couple hh and Ih components, it is convenient to 
use a single kind of oscillator states, e.g. hh states. This 
can be done by rewriting Eq. (|13|) in terms of the hh 
harmonic oscillator Hamiltonians: 



P-Q + Vc 



QD 



71 - 72 

7i +72 
71 + 2 72 
7i - 2 72 



Hhh,X 
Hhh. 
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71 +72 
272 



7i - 2 7^ 



Xx (x 2 + y 2 ) 
XzZ 2 , (14) 



where H hhj ± = T hh ,±_ - Xx (x 2 + y 2 )/2 and H h h,z = 
Thh.z~Xz {z 2 )/2. The resulting hole states are Luttinger 
spinors of the form: 

K >=E < O.K'«V' w *)l 3 /2 J J*>, (15) 

r,J z 

where Vj = 0, 1, 2 . . . is the quantum number of the ID 
hh harmonic oscillator along the j direction, r is the 
combined orbital quantum number, r — {v x ,v y ,v z ) and 
1 3/2 , J z ) the Bloch function. 

H p h in Eq. (Q} is the Hamiltonian of acoustic phonons, 
given by H ph = J] q A Sw «a »q A a qA, with u qX standing 
for the phonon energy spectrum of branch A (A = l,tl,t2 
for longitudinal and the two transversal phonon modes) 
and momentum q. We restrict to low phonon energies, 
where the linear dispersion regime applies, uj q \ = c\ q, 
with c\ as the phonon velocity. 

Hh-ph is the hole-phonon interaction, 



H h - 



h—ph 



e(p pz l+ H, 



dp • 



(16) 



where e is the hole charge, <j} pz the piezoelectric poten- 
tial and Hdp the deformation potential term. These 
are the two relevant scattering mechanisms at low 
temperatures^ The piezoelectric potential is given by>22 



b VZ = ^P* 



E 

Aq 



4-7T i 

e r q z 



'14 



(<?x e\ z 



1y £ zx 



1Z £ X y) 



(17) 



where e r is the relative dielectric constant, /114 the piezo- 
electric constant and £y the strain tensor component. On 
the other hand, the deformation potential term is given 
by the Bir-Pikus strain Hamiltonian: 



/ p X + q X 
-(a*)t 
( r A)t 



V 







fT - q* 


(r A )t 





p x - q x 





v 
s 



p X + q X ) 



(18) 
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where: 



for transversal tl phonons, and 



P = -a(4+4/+4)> 
b , , 



£ A ) 


(19) 


p' 2 


24), 


(20) 


r* 2 




(21) 
(22) 





r * = ^ 6 (e A - e A ' 
s —d(e zx — i£y Z ). 

Here a, 6 and d are the valence band deformation poten- 
tial constants. 

The components of the strain tensor are rewritten us- 
ing normal-modes coordinates, 29 

4 = ~\ E^frKq)^ +V A (q)^(q,r), (23) 



= 0, 



2 q A 



d TTti„\ ST ^ (lx~iq y ) 



s t2 = ~ U\q) ^ — F(q, r), (27) 



for transversal tl phonons. 



B. Relaxation rate 



where F(q,r) = {e~^ v a+ + e^ r a q ) and U x (q) = 
yJhj2pVujq\, with p and V standing for the crys- 
tal density and volume. ^a(q) are the phonon po- 
larization vectors: 7y/(q) = (q x , q y , q z ) /q, r] n (q) = 

(qxqz,q v qz,-ql) /qq± and r? t2 (q) = (q y ,-q x ,Q) /q±, 



with q± = y<7 2 +<7 2 - The piezoelectric potential now 
reads: 

6 r q 

q 



rtl = _47T/ll4 



E 



^(2 9z - g± ) 



6' 2 

/J 2 



4tt/ii4 t <?2 (g 2 - g 2 ) 

^ E — ^7 ^(q>r), (24) 



In turn, the deformation potential operators become: 
V l = iaU\q) ^ 9 F(q,r), 



3 = v 




;^ \ qz (q x -iq y ) 



s l = idU\q) Y, 
q 

for longitudinal phonons, 

P n = 0, 



F(q,r) 



(25) 



r 41 = -t 



y (^-b Sl^k^l ld «- * 9- 



2 gg_L 



^(q,r), 

s ii = „■ d L rf (;) E — = ^ ^ X " - — 

q q^q 



F(q,r), (26) 



We calculate the spin relaxation from an initial hole 

h\ 
fl 



state l^' 1 ), with energy Ef, to a final hole state 



with energy E'j. The relaxation rate is estimated with 
a Fermi golden rule. We consider zero temperature, so 
that there is no phonon absorption. After integrating 
over phonon degrees of freedom, the rate is given by: 



2tt 

T 



E 



^ h f\nt ph \^ 



5(AE fi + hc x q). 



(28) 

Here 'Hf^_ p f l is the hole-phonon interaction Hamiltonian, 
Eq. (|16|) , but for a fixed phonon branch A and momentum 
q, and AE fl = E h f - E£. It can be seen from Eqs. ((241 

[27f that all the terms of Hf l q _ p f l contain a factor which 
depends on q only and F(q, r), which depends on spatial 
coordinates as well. Thus, when expanded, the matrix 
Xq ^) takes the form: 



element ^>)\U h _ ph 



&)\< q - 



ph 



J' z ,J Z: r' ,r 



(29) 



where GV,r'(q) = e lqr \r) and Mj, j (q) gathers the 
q-dependent factor of the 'H Xq _ ph term coupling J z and 
J' z . G r y (q) is calculated analytically using iterative pro- 
cedures as described in Ref. l40l . The sum over q in 
Eq. (|28|) is then carried out using numerical integration. 
To this end, it is convenient to use spherical coordinates, 
as the modulus q is fixed by the resonance condition and 
we are left with a two-dimensional integral. 

Calculations are carried out for InAs QDs embedded 
in a GaAs matrix. When differences are expected, we 
also calculate GaAs QD embedded in an Al x Gai_a;As 
matrix. Table U summarizes the parameters we use. The 
parameters correspond to the QD material, except for the 
crystal density and velocity of sound, which correspond 
to the matrix material because we assume bulk phonons 
(for simplicity, for Al^Gai-^As we assume x — > and 
use GaAs phonon parameters). For the dielectric con- 
stant, an average value of e r = 12.9 is taken all over the 
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TABLE I: Parameters used in the numerical calculations for 
InAs (left column) and GaAs (right column) QDs. GaAs pa- 
rameters are used for the matrix in both cases. e~ , h + and 
ph stand for electron, hole and phonon. 



Parameter 


Symbol 


InAs 


GaAs 


R,cf. 


e mass (JTtg) 


7Tl e 


0.026 


0.067 


A 1 
41 


h Luttingcr param. 


71 


20 


6.98 


41 


Luttingcr param. 


72 


8.5 


2.06 


A 1 

4:1 


Luttingcr param. 


73 


9.2 


2.93 


41 


e~ deformation pot. (cV) 


a c 


-5.08 


-7.17 


41 


deformation pot. (cV) 


a 


1.0 


1.16 


41 


h~ >! ~ deformation pot. (cV) 


i. 
o 


-1.8 


-2.0 


41 


h' >r deformation pot. (cV) 


c 


-3.6 


-4.8 


41 


e~ BIA cocff. (cV-A 3 ) 


b 41 


27.18 


27.58 


28 


h + BIA cocff. (cV-A) 


c k 


-0.0112 


-0.0034 


28 


h+ BIA cocff. (cV-A 3 ) 


&41 


-50.18 


-81.93 


28 


h+ BIA cocff. (cV-A 3 ) 


b42 


1.26 


1.47 


28 


h+ BIA cocff. (cV-A 3 ) 


&51 


0.42 


0.49 


28 


h + BIA cocff. (cV-A 3 ) 


bs2 


-0.84 


-0.98 


28 


Longitudinal ph speed (m/s) 


Cl 


4720 


4720 


41 


Transversal ph speed (m/s) 


Ct 


3340 


3340 


41 


Crystal density (kg/m 3 ) 


P 


5310 


5310 


41 


Piezoelectric coeff. (V/cm) 


hl4 


3.5-10 6 


1.45-10 7 


49 



structure. The basis used to solve Hamiltonian ([T]) con- 
tains all the hh oscillator eigenstates with the quantum 
numbers v x ,v y < 13 and v z < 9. 

III. NUMERICAL RESULTS AND DISCUSSION 

We shall start this section by describing the depen- 
dence of the hole spin lifetime on the QD geometry and 
the spin splitting magnitude (Section MI A[) . The influ- 
ence of each parameter can be understood by analyzing 
the spin admixture mechanisms, as we show in Section 
IfflBl Next, in Section UTTCl we study the effect of the 
ground state changing from mainly hh character, which is 
the case addressed in most previous studies, to mainly Ih 
character. This transition is observed in QDs with large 
aspect ratio42~— Last, in Section fill D| we compare the 
role of deformation potential and piezoelectric potential 
interactions in determining the efficiency of hole-phonon 
coupling. 

A. Geometry and spin splitting dependence 

Solid lines in Figure [5] show the hole spin lifetime as 
a function of the QD geometry and the spin splitting 
energy of InAs QDs. For comparison, we also plot the 
electron spin lifetimes (dotted lines). The latter have 
been calculated using the same formalism as for holes 
but adapted for single-band conduction electrons^ Both 
electrons and holes are assumed to feel the same force 
constant, but we use the hh oscillator frequency, u) z h or 
uif 1 , to denote the confinement strength (for simplicity, 



we shall often drop the hh superscript). One can see im- 
mediately in the figure that the behavior of holes differs 
drastically from the well-known case of electrons. Below 
we summarize the influence of each parameter. 

Fig. [UJa) shows the spin lifetime dependence on the 
lateral confinement in QDs with strong vertical confine- 
ment. Tf increases monotonically with uj±. This is be- 
cause, as we approach the spherical confinement regime 
(lu z = u>j_), the Dresselhaus SOI of electrons is gradu- 
ally suppressed^ No such trend is however observed for 
holes, as Hbia does not cancel out even if the confine- 
ment is isotropic. As a matter of fact, T-f 1 shows an evi- 
dent non-monotonic behavior, with minimum at = 28 
meV and increasing away from it. It is worth noting that 
a previous study by Woods et al. predicted T/ 1 to de- 
crease with the QD diameter—, while a similar study by 
Lu and co-workers for somewhat larger QDs predicted 
the opposite trends Fig. Ufa) shows that both predic- 
tions are conciliable, corresponding to the right and left 
sides of the minimum, respectively. The origin of the 
different trends will be discussed in Section MI Bl 

Figs. [D^b-c) show the spin lifetime dependence on the 
vertical confinement in QDs with moderately strong and 
weak lateral confinement, respectively. The confinement 
strengths roughly correspond to self- assemble d 46 ' 47 and 
electrostatic^ QDs. As can be seen in Fig.[5Jb), electrons 
and holes have opposite behaviors. Tf now decreases 
with uj z because the structure is becoming flatter (less 
isotropic). Instead, the behavior of T/ 1 is similar to that 
observed for varying lateral confinement, with a shallow 
minimum at lo z — 14 meV. Previous studies have shown 
that T/ 1 increases with the vertical confinement. 30 This 
is consistent with the right side of the T± minimum in 
Fig. [Hb), but here we show that the opposite trend is 
also possible if the QD aspect ratio is large enough (left 
side of the minimum). 

Fig- IHc) illustrates the case in which the lateral con- 
finement is weak. The minimum of Ty is now shifted 
towards very small huj z values so that only the right side 
behavior is seen in the range under study. Interestingly, 
here T/ 1 shows a clear saturation with increasing vertical 
confinement (hw z > 40 meV), which has not been noticed 
before^ As we show in Section fill Bl this saturation re- 
flects the fact that Hbia has replaced hh-lh mixing as 
the main source of SOI. 

The inset in Fig. &b) compares T/ 1 in InAs and GaAs 
QDs with moderate lateral confinement. As can be ob- 
served, the spin lifetime in InAs QDs is longer than that 
in GaAs QDs when uj z > cj^, which is e.g. the case 
of self-assembled QDs. This is an unexpected result be- 
cause the valence band SOI of InAs is stronger than that 
of GaAsi 41 - The underlying reason is that the Dresselhaus 
SOI interaction is stronger for GaAs (see 641 in Table fij, 
and it provides a significant contribution to the total spin 
admixture. 

Figure [Hd) shows the the spin-flip lifetime of electrons 
and holes in a self-assembled-like QD as a function of the 
spin splitting energy. For electrons, Tf is largely deter- 



6 




A (meV) ftco z (meV) 

FIG. 2: Hole (solid line) and electron (dotted line) spin lifetime in InAs QDs, as a function of the lateral confinement (a), 
vertical confinement (b-c) and spin splitting energy (d). The inset in (b) compares Ti for InAs and GaAs. (a): hJl h = 50 
meV, A = 0.4 meV. (b): huf 1 = 20 meV, A = 0.4 meV. (c) hcuf 1 = 5 meV, A = 0.4 meV. (d): rvJf = 20 meV and Tilo^ = 50 
meV. 



mined by the efficiency of the carrier-phonon coupling. 
It is short when the phonon wavelength is of the same or- 
der as the carrier wavefunction extension, but it increases 
for large (small) A values because the phonon wavelength 
becomes too short (long), as q s» A/ he. The same hap- 
pens for holes (notice that G r y (q) in Eq. vanishes 
in the limits of large and small phonon wavevector). Yet, 
figure H{d) shows that T/ 1 is only sensitive to A for small 
splittings, but then it reaches a plateau where Ty ~ fis. 
The different behavior of holes and electrons is due to 
the different effective mass along the growth direction, 
m z . For hu> z h = 50 meV, the characteristic length of the 
oscillator states in the growth direction, l z — ^/hjm z lu z , 
is l z — 4.5 nm for electrons and l z h = 2.4 nm for holes, 
i.e. the hole confinement is stronger. As a result, larger 
values of A are required for T/ 1 to increase. 

Experiments with excitons in self-assembled InGaAs 
QDs have revealed a negligible influence of e-h inter- 
actions on T/\— Our results would be consistent with 
such a observation in the regime of large A. For self- 
assembled QDs, however, A < 0.5 meV. Thus, the insen- 
sitivity noticed in the experiment is inconsistent with the 
single-phonon processes we consider. This suggests that 
two-phonon processes dominate in these systems! 15 i 32 On 
the other hand, experiments with colloidal nanorods have 
shown a sizable increase of Tj when changing from type- 
I to type-II confinement, which modulates the e-h ex- 
change energy from typical colloidal structure values (few 
meV) to type-II system values (fractions of meV)jS We 
have run simulations for a nanorod-like geometry (dashed 
line in Fig.^d)) and find that the weak vertical confine- 



ment renders Ty sensitive to A for all the range under 
study, in agreement with the experiment. This indicates 
that the weak vertical confinement of rods renders single- 
phonon processes efficient. 

To close this section, we notice that previous theoret- 
ical studies have predicted that hole spin lifetimes can 
exceed those of electrons in very flat QDs. 31 Fig.[5]shows 
that this could actually be achieved in gated structures, 
where lateral confinement is very weak (see the crossing 
between and Tj in panel (c)). However, for typical 
self- assembled InAs QDs, Tf is about one order of magni- 
tude longer than if (see panel (b) for large fuo z )£2. This 
is consistent with experimental measurements by Heiss 
et al^S 



B. Mechanism of spin admixture 

Spin admixture is a requirement for spin-flip transi- 
tions to take placed In this section, we compare the 
spin admixture resulting from all the SOI terms affecting 
the hole ground state. As we shall see, once the domi- 
nant SOI mechanism is determined, one can rationalize 
the geometry dependence of T/ 1 described in the previous 
section. For convenience of the analysis, in what follows 
we plot and discuss relaxation rates, 1/T/ 1 . 

Figure [3Ja) shows the relaxation rate for the InAs 
QDs of Fig. [2ja), but now obtained by calculating hole 
states with the diagonal terms of Hl plus different SOI 
terms: off-diagonal Hl terms (hereafter hh-lh coupling), 
full Dresselhaus Hamiltonian (Hbia), hh-lh coupling plus 
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linear-in-fc term (Hl + Hc k ) and hh-lh coupling plus the 
dominant cubic- in-fc Dresselhaus term (Hl +^6 4 i)- The 
total rate, corresponding to Hl + Hbia, is also shown 
(thick black line). One can see that Hl (red solid line) is 
more important than Hbia (blue dashed line) for large 
lj±_ values. However, as the lateral confinement is weak- 
ened and the system becomes flatter, Hbia gains weight. 
For self-assembled QDs (huij_ ss 10 — 25 meV), Hbia 
is already comparable to Hl and it becomes dominant 
for weakly confined (e.g. gated) QDs. Fig. [3] also re- 
veals that the linear-in-fc BIA term (gray dot-dashed 
line) is but a secondary mechanism, which barely en- 
hances the relaxation rate coming from Hl- This is in- 
spite of the fact that it is a source of direct admixture 
between hh states with opposite spin, with no participa- 
tion of Ih states. For this reason, it had been proposed 
as the dominant SOI term in flat QDs with weak lateral 
confinement.^ Instead, most of the Hbia contribution 
comes from the cubic-in-fc 641 term (see green dotted 
line). This term relies on intermediate Ih states in or- 
der to mix hh J z = +3/2 and J z = —3/2 components 
(see Eq. (|A3[) in the Appendix), which implies that a si- 
multaneous description of hh and Ih states is necessary 
for realistic modeling. 



,nn , ^ 
l x (nm) 

15 10 8 7 6 5 




10 20 30 40 50 60 70 80 90 
riC0 ± (meV) 



FIG. 3: Hole spin relaxation (a) and weight of the minor 
hh component (b) as a function of the lateral confinement. 
Different SOI terms are considered. The system is the same 
as that of Fig. [2J a). 

By comparing the total relaxation rate coming from 
Hl + Hbia with that coming from Hl and Hbia, it is 
clear that the total rate is not just the sum of the two 
independent mechanisms. For example, at fkoj_ = 30 
meV, adding Hbia to Hl enhances 1 /Xf about an order 



of magnitude, even though the contribution coming from 
Hbia alone is about 100 times smaller than that coming 
from Hl- Thus, simultaneously accounting for both SOI 
terms is also required for quantitative estimates. 

For a qualitative understanding of the geometry de- 
pendence of 1/T-l , we rewrite the hole states, Eq. (|15p . 
as |#* >= E^c£|<^>|3/2,J Z >, where |^> = 
J2 r c ™-r \ r ) i s the envelope function associated with the 
periodic Bloch function |3/2, J z ). If we restrict to the 
diagonal components of Hh- P h, the matrix element de- 
termining the relaxation rate becomes: 

(*f\< g - ph l*?> « E ( c i ^jM- p h Itfr.)- (30) 

When the QD is oblated or spherical, the low energy 
states are essentially hh states. Thus, the initial (fi- 
nal) state of the spin transition is mostly a pure spin up 
(spin down) hh state. Considering that \cy 2 \ 3> I c 3/2 1 

d c -3/2l l c -3/2l)' one can obtain an approximate ex- 
pression: 




2 



(31) 

In other words, the relaxation rate is proportional to the 
spin admixture of the ground state through the squared 
coefficient of the minor hh component (here spin up, 
J z = 3/2), and proportional to the efficiency of the hole- 
phonon coupling through the rightmost matrix element. 

The geometry dependence of 1/T{ 1 for a given SOI 
mechanism simply reflects the spin admixture variation. 
This can be seen in Fig. [3fb) , which shows that, for Hl 
and Hbia, I £3/2 1 2 has the same qualitative dependence 
on the geometry as the corresponding 1/T/ 1 values in 
Fig. EJa). This allows us to interpretate the observed 
maximum as a function of uj±. For Hl, the spin 
admixture takes place necessarily through Ih states, 
see Eq. The weight of the minor hh component is 
related to the hh-lh coupling strength (through R and 
S operators in Hl) and to the energy splitting between 
the hh and intermediate Ih states. The hh-lh coupling 
terms scale up with the confinement strength (note R 
and S operators dependence on k). In turn, the energy 
splitting, is largely influenced by the different masses 
of hh and Ih. It is maximum for flat structures and 
goes to zero for slightly prolate structures. 44 ' 45 On the 
left side of Fig. [3l uj± <C cj z , so Ih states are far in 
energy. Changing u>± has a weak effect on the energy 
splitting, which is mainly set by the lateral confinement. 
Then, the main effect of incresing oj± is to enhance the 
hh-lh coupling strength, resulting in larger |c^ 2 | values. 
Around huj_ rs 30 meV, however, the lateral confinement 
energy becomes comparable to the vertical one. The Ih 
states get close in energy and start participating in the 
admixture themselves, thus reducing the spin admixture 
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of the minor hh component. 



i hh , n 
1 z (nm) 




40 60 
hco z (meV) 



FIG. 4: Hole spin relaxation as a function of the vertical con- 
finement. Different SOI terms are considered, (a): strong 
lateral confinement, hu)± — 20 meV. (b): weak lateral con- 
finement, fajj± = 5 meV. The inset in (a) is the corresponding 
result for GaAs. The system is the same as that of Figs. [2jb) 
and (c). 



Next, we analyze the mechanisms involved in the spin 
relaxation with varying vertical confinement. Figure 0] 
shows 1 /Ti for the same systems as in Fig. [2jb-c) , but 
calculating the hole states with the diagonal terms of Hl 
plus hh-lh coupling (red solid line) or full Hbia Hamilto- 
nian (blue dashed line). Panel (a) corresponds to strong 
lateral confinement. The total rate has a maximum at 
fku z — 14 meV, whose origin is analogous to that de- 
scribed above for varying lateral confinement. Beyond 
the maximum, the total rate {Hh + Hbia) decreases 
with the vertical confinement strength, in agreement with 
previous studies^ Yet, we also find that the decrease 
eventually saturates. This is evident for InAs QDs with 
weak lateral confinement, as shown in Fig.UJb), or GaAs 
QDs with strong lateral confinement, as in Fig. |4ja) in- 
set. The origin of this saturation is the contribution of 
Hbia, which provides a lower bound to l/T^ . In partic- 
ular, i?b 41 introduces off-diagonal coupling terms which 



are quadratic in k z (see operator L41 in the Appendix), 
instead of the linear k z terms of Hl (see S operator in 
Eq. (@}). Since the uncoupled hh and Ih energies are also 
quadratic in k z , a perturbational analysis easily shows 
that the two contributions compensate each other. For 
strong uj z , when lateral confinement is negligible, the can- 
cellation is exact and the relaxation rate does not depend 
on the vertical confinement. 




FIG. 5: Hole spin relaxation as a function of the spin splitting 
energy in a QD with huu±_ = 20 meV and hio z = 50 meV. Same 
legend as in Fig. [3] is used. The inset compares the energy 
splitting between the Kramer doublet for hh-lh coupling and 
Dresselhaus SOI. 

The magnitude of the spin splitting also influences the 
dominant mechanism of spin admixture. This is illus- 
trated in Fig. 03 which considers a self-assembled InAs 
QD with varying spin splitting energy. For large A, Hl 
has a dominant contribution to the relaxation rate, but 
Hbia becomes equally important for small enough A. 
The relative enhancement of the role of Hbia originates 
in its zero- field spin splitting, which leads to larger energy 
difference between the Kramers pair (AEfi) than with A 
alone, as shown in the figure inset. When A — > and the 
phonon wavelength increases beyond the QD size, the 
extra energy coming from the zero-field spin splitting of 
Hbia provides a significant contribution to preserve the 
hole-phonon coupling efficiency. 



C. Light hole spin relaxation 

When the aspect ratio increases and the QD shape 
becomes prolate, the hole ground state evolves from the 
eminent hh character discussed so far to Ih character, as 
noted e.g. in colloidal nanorods42r— Here we investigate 
how the change of the ground state affects the relaxation 
between the two highest hole states. Fig. [D^a) shows the 
hole energy levels in a QD with haj± = 40 meV as a 
function of hio z . In the limit of strong and weak vertical 
confinement, the two highest states are essentially hh and 
Ih doublets, J z = ±3/2 and J z = ±1/2 respectively. In 
the intermediate regime, hu> z =9 — 17 meV, the two 
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doublets cross. This gives rise to pronounced changes in 
the relaxation rate, as shown in Fig.[SJb). 
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FIG. 6: Hole energy levels (a) and spin relaxation rate (b) 
around the hh-lh crossing region. fyjJj_ = 40 meV and A = 2 
meV. In (a), solid and dashed lines are used for states with 
dominant Ih and hh character, respectively. The dotted line 
gives the energy splitting. Shades are used to identify the 
regions with different kind of states involved in the transition. 

The changes can be understood as follows. In region I 
of the figure, the two highest states are the hh doublet. 
The relaxation is between states with opposite spin and 
roughly constant energy splitting (see AEfi, dotted line 
in Fig. |6](a)). At ftu z = 17 meV, when we enter region 
II, the excited hh state crosses with the first Ih state. 
Now the relaxation is between a Ih (J z — —1/2) and a 
hh (J z = —3/2). Since Ih have mixed spin up and spin 
down projections, there is no need for spin flip and the 
transition is much faster. This explains the abrupt in- 
crease of l/T^ . However, the energy splitting between 
the hh and Ih becomes smaller with decreasing uj z be- 
cause of their different masses. As a result, at Hu> z = 11 
meV, the Ih replaces the hh as the ground state. Near 
the degeneracy point, haj z = 11.2 meV, AEfi is so small 
that hole-phonon coupling becomes inefficient and the re- 
laxation is strongly suppressed, but it increases again in 
region III for the same reasons as in region II. Finally, at 
hoj z = 9 meV, the excited Ih state crosses with the high- 
est hh state. We thus enter region IV, where the tran- 
sition takes place between two Ih states with orthogonal 
Bloch functions, 1 3/2, ±1/2). J z admixture mechanisms 
are necessary and the relaxation becomes slow. As a mat- 
ter of fact, the spin relaxation timescale for transitions 
between Ih states is similar to that between hh states, in 
spite of the fact that their Bloch functions contain mixed 



spins. 



D. Mechanism of hole-phonon coupling 

Electron-acoustic phonon coupling in QDs is known to 
occur mainly through deformation potential interaction 
when the energy splitting is large (AEfi ^ 0.1 meV) 
and through piezoelectric potential when it is small^ In 
principle, for holes the situation may differ because the 
deformation potential Hamiltonian, Eq. (|18[) . is formally 
different from that of electrons. To investigate this point, 
in this section we compare the role of the two kinds of 
carrier-phonon coupling mechanism for holes subject to 
varying effective magnetic fields. 
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FIG. 7: Hole spin relaxation as a function of the spin splitting 
energy, (a): transition between hh states in a QD with hw±_ = 
20 and fkj z — 50 meV. (b): transition between Ih states in a 
QD with hw± — 40 and hw z — 5 meV. 

Figure Efa) shows the spin relaxation rate in an oblate 
(quasi-2D) QD, where the highest states are hh, while 
Fig. EJb) shows that in a prolate (quasi-lD) QD, where 
the highest states are Ih. For the spin transition within 
the hh doublet, panel (a), the behavior is analogous to 
that of electrons. Deformation potential interaction (dot- 
ted line) provides the largest contribution to 1 /T± except 
for very small A, when piezoelectric interaction (dashed- 
dotted line) takes over. This is because all the terms of 
Hdp are proportional to the phonon momentum J~q, see 
Eqs. (|25H27|) . while the piezoelectric potential is propor- 
tional to 1/y/q, see Eq. (f24|) . With decreasing A both 
mechanisms become inefficient, because for long phonon 
wavelength e lqr ps 1. Then, in Eq. (|29|) . the matrix el- 
ement G r y(q) w (r'\r), i.e. it tends to the overlap be- 
tween the envelope components of the initial and final 
states. For hh, these components have different symme- 
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tries, so the coupling vanishes. For example, in axially 
symmetric structures, the J z = +3/2 component of the 
initial state has a well defined azimuthal angular momen- 
tum m z = 0, which couples through the s A operator of 
Hdp with the J z = +1/2 component of the final state, for 
which m z = — 

The situation for Ih is quite different. As shown in 
Fig. EJb), now deformation potential interaction is the 
dominant relaxation channel even for small A. The un- 
derlying reason is that, in contrast to the hh case, the 
off-diagonal terms of Hd P couple envelope components 
with the same symmetry. For example, the J z = +3/2 
and J z = +1/2 components of the initial and final state 
now have both m z — — 1, and hence are not orthogonal. 
As a consequence, G r . r '(q) does not vanish when q — » 0. 



Appendix A: Dresselhaus Hamiltonian for Holes 

In this appendix we write the explicit matrix forms of 
the different Hbia terms in cartesian coordinates. Sep- 
arating the different invariants in Eq. @ we obtain: 



Hbia — Hc k + Hb iX + Hb B1 + Hb S2 



where: 



H, 



Cfc 



= c. 
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-k z 
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(Al) 



(A2) 



IV. SUMMARY 



with k± — k x ± i k y . 



We have investigated hole spin relaxation in InAs 
and GaAs QDs using 3D 4-band k-p Hamiltonians. 
We have shown that the hole spin lifetime has a non- 
monotonic dependence on the lateral and vertical con- 
finement strength. This is due to the interplay between 
the energy splitting of hh and Ih, which is set by their dif- 
ferent masses, and the anisotropic hh-lh coupling terms. 
The resulting geometry dependence of hole spin relax- 
ation is qualitatively different from that of electrons. 

hh-lh coupling and Dresselhaus SOI have been found 
to have comparable contributions to the spin admix- 
ture of hole states in self-assembled QDs, with the for- 
mer becoming dominant for prolate structures, such as 
nanorods, and the latter for strongly oblate ones, such as 
gated QDs. The cubic-in-fc Dresselhaus term leads to an 
upper bound of T[ l with increasing vertical confinement, 
which is missed when considering hh-lh coupling only. 

The spin relaxation between hh and Ih is very fast be- 
cause Ih states have strong spin admixture. Instead, spin 
transitions between Ih states with different Bloch angular 
momentum J z take place in similar time scales as those 
between hh states with orthogonal spin. The main differ- 
ence is that, for Ih, deformation potential is the principal 
hole-phonon interaction mechanism even for small energy 
splittings. 
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where K + = K x + i K y , K_ =K x -i K y , K x = k x (k? y 
^z)> ~ k y (k x + k z ), and K z = k z (k x + k y ). 
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